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Abstract: We study a three dimensional continuous model of gravitating mat- 
ter rotating at constant angular velocity. In the rotating reference frame, by 
a finite dimensional reduction, we prove the existence of non radial stationary 
solutions whose supports are made of an arbitrarily large number of disjoint 
compact sets, in the low angular velocity and large scale limit. At first order, 
the solutions behave like point particles, thus making the link with the relative 
equilibria in A^-body dynamics. 
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1. Introduction and statement of the main results 

We consider the Vlasov-Poisson system 

dtf + v- V,/ - V,(/) • V„/ = 



1 



47r 



* P , P / f dv 



(1) 



R3 



which models the dynamics of a cloud of particles moving under the action of 
a mean field gravitational potential solving the Poisson equation: A(f> — p. 
Kinetic models like system ([T]) are typically used to describe gaseous stars or 
globular clusters. Here / = f{t,x,v) is the so-called distribution function, a 
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nonnegativc function in L°°(R,i^(K'^ x M'^)) depending on time t G M., position 
X gM.^ and velocity v £ M^^, which represents a density of particles in the phase 
space, M.^ x R"^. The function p is the spatial density function and depends only 
on t and x. The total mass is conserved and hence 



does not depend on t. 

The first equation in ([T]) is the Vlasov equation, also known as the collisionless 
Boltzmann equation in the astrophysical literature; see [S]- It is obtained by 
writing that the mass is transported by the flow of Newton's equations, when 
the gravitational field is computed as a mean field potential. Reciprocally, the 
dynamics of discrete particle systems can be formally recovered by considering 
empirical distributions, namely measure valued solutions made of a sum of Dirac 
masses, and neglecting the self-consistent gravitational terms associated to the 
interaction of each Dirac mass with itself. 

It is also possible to relate ([T]) with discrete systems as follows. Consider 
the case of N gaseous spheres, far away one to each other, in such a way that 
they weakly interact through gravitation. In terms of system ([1]) , such a solution 
should be represented by a distribution function /, whose space density p is 
compactly supported, with several nearly spherical components. At large scale, 
the location of these spheres is governed at leading order by the TV-body gravi- 
tational problem. 

The purpose of this paper is to unveil this link by constructing a special class 
of solutions: we will build time-periodic, non radially symmetric solutions, which 
generalize to kinetic equations the notion of relative equilibria for the discrete 
A^-body problem. Such solutions have a planar solid motion of rotation around 
an axis which contains the center of gravity of the system, so that the centrifugal 
force counter-balances the attraction due to gravitation. Let us give some details. 

Consider N point particles with masses ruj , located at points Xj (t) £ M"^ and 
assume that their dynamics is governed by Newton's gravitational equations 



Let us write a; e R-^ as a; = S x R « C x M where, using complex 

notations, x' = {x^,x'^) « x^+ix"^ and rewrite system ^ in coordinates relative 
to a reference frame rotating at a constant velocity a; > around the x'^-axis. 
This amounts to carry out the change of variables 




p[t, x) dx = M 




N 




z' =^ z^ + i z' 



,2 



In terms of the coordinates {z',z^), system ^ then reads 




N 



dt^ 



E 



TOfc Zk - Zj 



A-n \zk - Zj\^ 



j = l,...iV. (3) 
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Wc consider solutions which are stationary in the rotating frame, namely con- 
stant solutions (zi, . . . zn) of system ([3]). Clearly all Zj's have their third com- 
ponent with the same value, which we assume zero. Hence, we have that 

^fc = (6,0), &eC, 

where the ^fc's are constants and satisfy the system of equations 

N 

+ ^'?, , j- = l,...iV. (4) 



In the original reference frame, the solution of ^ obeys to a rigid motion 
of rotation around the center of mass, with constant angular velocity uj. This 
solution is known as a relative equilibrium, thus taking the form 

x^{t)^{e'^'^,,0) , eC, j = l,...iV. 

System ^ has a variational formulation. In fact a vector (d, . . .^at) solves ^ 
if and only if it is a critical point of the function 



^ X 1 \ - Too mfc UJ ir-^ _ 



Here to denotes {mj)^^^. A further simplification is achieved by considering the 
scaling 

V^(Ci,...eiv) = 6.2/3 v,„(Ci,...Civ) (5) 

where 

AT N 



^ N ^ N 



This function has in general many critical points, which are all relative equilibria. 
For instance, Vm clearly has a global minimum point. 

Our aim is to construct solutions of gravitational models in continuum me- 
chanics based on the theory of relative equilibria. We have the following result. 

Theorem 1. Given masses mj, j = 1, . . . N , and any sufficiently small lu > 0, 
there exists a solution fuj{t, x, v) of equation ([T} which is ^ -periodic in time and 
whose spatial density takes the form 

N 

p(t,x) := / dv ^Y.Pi^'' - W) + • 

Here o(l) means that the remainder term uniformly converges to as lu ^ 0+ 
and identically vanishes away from UjLiBii{x'j{t)), for some R > 0, independent 
of CO. The functions Pj{y) are non-negative, radially symmetric, non-increasing, 
compactly supported functions, independent of uj, with J^g^s Pj{y) dy = nij and 
the points x'j{t) are such that 

x^{t)=Lj-^/^{e'^'a,0) , aeC, j^l,...N 
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and 



lim V™ (Cr , . . . = min V„ , lim VV,„ (C^ ,...(%)= 



The solution of Theorem [T] has a spatial density which is nearly spherically 
symmetric on each component of its support and these ball-like components 
rotate at constant, very small, angular velocity around the .T'^-axis. The radii of 
these balls are very small compared with their distance to the axis. Wc shall call 
such a solution a relative equilibrium of ([1]), by extension of the discrete notion. 
The construction provides much more accurate informations on the solution. 
In particular, the building blocks pj are obtained as minimizers of an explicit 
reduced free energy functional, under suitable mass constraints. 

It is also natural to consider other discrete relative equilibria, namely critical 
points of the energy Vm that may or may not be globally minimizing, and ask 
whether associated relative equilibria of system ([T]) exist. There are plenty of 
relative equilibria of the A^-body problem. For instance, if all masses mj are 
equal to some > 0, a critical point is found by locating the Q^s at the 
vertices of a regular polygon: 



r e 



2iTvij-l)/N 



j = l,...N, (6) 



where r is such that 



d_ 
dr 



ajv 1 2 

47r r 2 . 



with ajv := 



V2 ^ Vl-cos(27rj7iV) 



I.e. 



1 /3 

r = {ai\! m^/^in)) . This configuration is called the Lagrange solution, 
see j35j . The counterpart in terms of continuum mechanics goes as follows. 

Theorem 2. Let (Ci, . . . Cjv) be a regular polygon, namely with C,j given by 
and assume that all masses are equal. Then there exists a solution exactly as 
in Theorem[l\ but with \im^^o_^ {(f , . . . C,'^) = (Ci, • • ■ Cn)- 

Further examples of relative equilibria in the A^-body problem can be obtained 
for instance by setting — 1 points particles of the same mass at the vertices 
of a regular polygon centered at the origin, then adding one more point particle 
at the center (not necessarily with the same mass), and finally adjusting the 
radius. Another family of solutions, known as the Euler-Moulton solutions is 
constituted by arrays of aligned points. 

Critical points of the functional Vm are always degenerate because of their 
invariance under rotations: for any a G M we have 

V™(Ci,...Civ)=V„,(e*"Ci,---e'"Cjv). 

Let C = (Ci 7 • ■ • Ov) be a critical point of Vm with C,t ^ 0. After a uniquely defined 
rotation, we may assume that Q2 — 0. Moreover, we have a critical point of the 
function of 2A^ — 1 real variables 

V„(Ci,...Ctt,...Civ) := V,„(Ci,...(Cn,0),...Civ) ■ 
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Wc shall say that a critical point of V„, is non- degenerate up to rotations if the 
matrix Z)^V„i(^i, . . . (^i, . . . C,n) is non-singular. This property is clearly indepen- 
dent of the choice of ^. 

Palmore in [5nil5Tl[52 11331151] has obtained classification results for the rela- 
tive equilibria. In particular, it turns out that for almost every choice of masses 
mj, all critical points of the functional Vm are non-degenerate up to rotations. 
Moreover, in such a case there exist at least [2^~^{N — 2) -I- 1] (TV — 2) ! such 
distinct critical points. Many other results on relative equilibria arc available in 
the literature. We have collected some of them in [Appendix A| with a list of rel- 
evant references. These results have a counterpart in terms of relative equilibria 
of system ([T|). 

Theorem 3. Let (Ci, • • -Cjv) be a non-degenerate critical point of Vm up to ro- 
tations. Then there exists a solution f^ as in Theorem]^ which satisfies, as in 
Theorem^ lim^^o JCi" , ■ • • C5^) = (Ci,---Civ). 

This paper is organized as follows. In the next section, we explain how the 
search for relative equilibria for the Vlasov-Poisson system can be reduced to 
the study of critical points of a functional acting on the gravitational potential. 
The construction of these critical points is detailed in Section [31 Sections H] 
and [5] are respectively devoted to the linearization of the problem around a 
superposition of solutions of the problem with zero angular velocity, and to the 
existence of a solution of a nonlinear problem with appropriate orthogonality 
constraints depending on parameters (Cj)jLi related to the location of the N 
components of the support of the spatial density. Solving the original problem 
amounts to make all corresponding Lagrange multipliers equal to zero, which is 
equivalent to find a critical point of a function depending on this is the 

variational reduction described in Section [51 The proof of Theorems [1] [51 and [31 
is given in Section [71 while known results on relative equilibria for the A'^-body, 
discrete problem are summarized in [Appendix A 



2. The setup 

Guided by the representation ^ of the A^-body problem in a rotating frame, 
we change variables in equation ([1]), replacing x = (x' ,x^) and v = {v' ,v'^) 
respectively by 

(e*'^*x',a;3) and (i a;' -f e"^* w', w^) . 
Written in these new coordinates, Problem ((TJ becomes 



U ^ - .\ I * P , P'^ f f '^^ ■ 
47r • ./b3 



(7) 



The last two terms in the equation take into account the centrifugal and Coriolis 
force effects. System ([7]) can be regarded as the continuous version of problem (jSj. 
Accordingly, a relative equilibrium of System ([T]) will simply correspond to a 
stationary state of ([7]). 
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Such stationary solutions of ([7]) can be found by considering for instance 
critical points of the free energy functional 

II P{f)dxdv + \ II {\v\'-u:^\x'\')fdxdv-\ I iVC/pdx 

for some arbitrary convex function /?, under the mass constraint 

f dxdv^ M . 

A typical example of such a function is 



// 



/?(/) - \ «r' /' (8) 

for some q S (1, oo) and some positive constant k^, to be fixed later. The corre- 
sponding solution is known as the solution of the polytropic gas model, see [SJUl 

[Mug. 

When dealing with stationary solutions, it is not very difficult to rewrite the 
problem in terms of the potential. A critical point of under the mass constraint 
//r3xR3 f dxdv = M is indeed given in terms of U by 

fix,v)=^{\ + ^\v\' + Uix)-^u;'\xf) (9) 

where 7 is, up to a sign, an appropriate generalized inverse of /3'. In case (|H1), 

7(s) = {—s)^^'^~^\ where s+ = (s + \s\)/2 denotes the positive part of s. 
The parameter A stands for the Lagrange multiplier associated to the mass con- 
straint, at least if / has a single connected component. At this point, one should 
mention that the analysis is not exactly as simple as written above. Identity ([9]) 
indeed holds only component by component of the support of the solution, if this 
support has more than one connected component, and the Lagrange multipliers 
have to be defined for each component. The fact that 



|x|->oo 47r|a;| 

is dominated by —51^^ as \x'\ — >■ 00 is also a serious cause of trouble, which 
clearly discards the possibility that the free energy functional can be bounded 
from below if a; ^ 0. This issue has been studied in in the case of the so-called 
flat systems. 

Finding a stationary solution in the rotating frame amounts to solving a 
non-linear Poisson equation, namely 

AU = g {X + U{x) - ^Lu^\x'\^) ifxGsupp(p) (10) 

and AU = otherwise, where g is defined by 
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Hence, the problem can also be reduced to look for a critical point of the func- 
tional 

J[U]-=\ I \^U\^dx+( G{\ + U{x)- ^uj^\x'\^) dx- [ Xpdx, 

where A = A[a;,C/] is now a functional which is constant with respect to x, 
with value A^, on each connected component Ki of the support of p{x) = 
g [X[x, U] + U{x) — ^ up' and implicitly determined by the condition 

g (Ai + U{x) — \uP' |a;'|^) dx = rrii . 

Ki 

By G, we denote a primitive of g and the total mass is M ~ X^iLi Hence we 
can rewrite as 

1 /■ ^ r /■ 

J[U] = - |Vf/pdx + y G{X^ + U{x)- ^uj^\x'\^) dx-m^X, 

We may also observe that critical points of F correspond to critical points of 
the reduced free energy functional 



Q[p] 



I (h{p)-\u^\x'\^p) dx^\( iVC/pdx 



acting on the spatial densities if h{p) = g~^{—s) ds. Also notice that, using 
the same function 7 as in ([9|), to each distribution function /, we can associate 
a local equilibrium, or local Gibbs state, 

Gf{x,v)=j{pipix)) + ^\v\') 

where p. is such that g{p) = p. This identity defines p — p{p) ~ 9^^{p) a-s a 
function of p. Furthermore, by convexity, it follows that J^[f \ > J''[Gf] = G[p] if 
p{x) = f{x, v) dv, with equality if / is a local Gibbs state. See for more 
details. 

Summarizing, the heuristics are now as follows. The various components Ki 
of the support of the spatial density p of a critical point are assumed to be far 
away from each other so that the dynamics of their center of mass is described 
by the A^-body point particles system, at first order. On each component Ki, the 
solution is a perturbation of an isolated minimizer of the free energy functional J- 
( without angular rotation) under the constraint that the mass is equal to to^. In 
the spatial density picture, on Ki, the solution is a perturbation of a minimizer 
of the reduced free energy functional Q. 

To further simplify the presentation of our results, we shall focus on the model 
of polytropic gases corresponding to In such a case, with p := + |, 5 is 
given by 

gii^) = 

if the constant Kq is fixed so that Kg = 4^72 /„+°° (i+t) 9^ dt = (Zir)^ ttttzxTi- 
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Free energy functionals have been very much studied over the last years, 
not only to characterize special stationary states, but also because they provide 
a framework to deal with orbital stability, which is a fundamental issue in the 
mechanics of gravitation. The use of a free energy functional, whose entropy part, 
//r3xr3 (^(Z) dxdv is sometimes also called the Casimir energy functional, goes 
back to the work of V.I. Arnold (see [TJ[HI3S]). The variational characterization 
of special stationary solutions and their orbital stability have been studied by 
Y. Guo and G. Rein in a series of papers T^.TIlIIHlIinilMllSZllMlISi and by many 
other authors, see for instance [i[IU}, 22 .23. 2-Cl[g5 1 |in i llI l lil] . 

The main drawback of such approaches is that stationary solutions which 
are characterized by these techniques are in some sense trivial: radial, with a 
single simply connected component support. Here we use a different approach 
to construct the solutions, which goes back to [13] in the context of Schrodinger 
equations. We are not aware of attempts to use dimensional reduction coupled 
to power-law non-linearities and Poisson force fields except in the similar case 
of a nonlinear Schrodinger equation with power law nonlinearity and repulsive 
Coulomb forces (see [5]), or in the case of an attractive Hartree-Fock model (see 
|21j). Technically, our results turn out to be closely related to the ones in [Sl[7]. 

Compared to previous results on gravitational systems, the main interest of 
our approach is to provide a much richer set of solutions, which is definitely of 
interest in astrophysics for describing complex patterns like binary gaseous stars 
or even more complex objects. The need of such an improvement was pointed 
for instance in [20j . An earlier attempt in this direction has been done in the 
framework of Wasserstein's distance and mass transport theory in |27] . The point 
of this paper is that we can take advantage of the knowledge of special solutions 
of the iV-body problem to produce solutions of the corresponding problem in 
continuum mechanics, which are still reminiscent of the discrete system. 



3. Construction of relative equilibria 

3.1. Some notations. We denote by a; = (x', x^) e x M a generic point in M^. 
We may reformulate Problem pU)) in terms of the potential u ^ — U as fol- 
lows. Given N positive numbers Ai,...AAr and a small positive parameter cj, 
we consider the problem of finding N non-empty, compact, disjoint, connected 
subsets Ki of K'^, i = 1, 2 . . . A^, and a positive solution u of the problem 

TV 
1=1 

lim u{x) = , (12) 

|a:|— ^oo 

where Xi denotes the characteristic function of Ki. We define the mass and the 
center of mass associated to each component by 

:= / Pi dx and x'^ := — / x pi dx . 

We shall find a solution of (|11D as a critical point of J [u] = u] in case of ([5]), 
with the notations of Section [21 that is when G{—s) = for any s <E R. 
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The functional J reads 

JM = -/ \Vu\^dx-—-Y, {u^K + \Lo^\x'\''Y^\ax. (13) 

Heuristically, our method goes as follows. We first consider the so-called basic 
cell problem: we characterize the solution with a single component support, when 
u) = Q and then build an ansatz by considering approximate solutions made of 
the superposition of basic cell solutions located close to relative equilibrium 
points, when they are far apart from each other. This can be done using the 
scaling invariance, in the low angular velocity limit w — > 0+. The proof of our 
main results will be given in Sections |4H2l It relies on a dimensional reduction 
of the variational problem: we shall prove that for a well chosen u, J [u] = 

J2iLi^i~^^* ~ Kn(?i'---^w) + 0(1) for some constant e*, up to o(l) terms, 
which are uniform in a; > 0, small. Hence finding a critical point for J will be 
reduced to look for a critical point of \>^^^ as a function of (^1, . . . Cat). 



3.2. The basic cell problem. Let us consider the following problem 

- Aw = {w- If^ in . (14) 



Lemma 1. Under the condition ^ini^aii^oo w{x) = 0, Equation (jl4p has a unique 
solution, up to translations, which is positive and radially symmetric. 

Proof. Since p is subcritical, it is well known that the problem 

-AZ = ZP inBi(O) 

with homogeneous Dirichlet boundary conditions, Z = 0, on dBi{Q), has a 
unique positive solution, which is also radially symmetric (see jl5)). For any 
i? > 0, the function Zr{x) := /^^^/(p-i) z{x/R) is the unique radial, positive 
solution of 

-AZr = ZP inSfl(O) 

with homogeneous Dirichlet boundary conditions on dBfi{0). According to |14l 
115] . any positive solution of (jl4p is radially symmetric, up to translations. Find- 
ing such a solution w of (|14p is equivalent to finding numbers R > and 
> such that the function, defined by pieces as w = Zj^ + 1 in Bf> and 
w{x) = m^,/{ATT\x\) for any x g M'^ such that |.t| > R, is of class C^. These 
numbers are therefore uniquely determined by 

w{R-) = \ = ^= w{R+) , w'{R-) = i?- Z'(l) = --^ = w'{R+) , 



which uniquely determines the solution of ((E]). 



□ 
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Now let us consider the slightly more general problem 

-Aw^ = {w^ - X)\ in 

with \\m\^\^^w^{x) 0. For any A > 0, it is straightforward to check that it 
has a unique radial solution given by 

Let us observe, for later reference, that 

/ (u;^ - A)^ dx = A(3-P)/2 [ {w- 1)1 dx =: X^^-p)/^ . (15) 

Moreover, is given by 

y,>^(x) = a(3-p)/2 V .t G such that \x\ > RX^^p-^^'^ . 

47r|x| 

3.3. The ansatz. We consider now a first approximation of a solution of piD - 
(|12p , built as a superposition of the radially symmetric functions w^^ translated 
to points ^i, i = 1, . . . in X {0}, far away from each other: 

N 

Recall that we are given the masses mi, . . .rriM- We choose, according to for- 
mula (|15p. the positive numbers Aj so that 

/ {wi — Xi)^ dx = rrii for alH = 1, . . . A'^ . 

By ^ we denote the array (^i, ^2, ■ • ■ £,n)- 

We shall assume in what follows that the points £,i are such that for a large, 
fixed /i > 0, and all small a; > we have 

le,|<M^-'/^ ie.-oi>M-'c.-2/'. (16) 

Equivalently, 

IC«I<M, \Q^Q\>fi-' where Q := uj-^/'' . (17) 
We look for a solution of (fTTj) of the form 

u = + 

for a convenient choice of the points ^i, where (j) is globally uniformly small when 
compared with W^. For this purpose, we consider a fixed number R > 1 such 

that 

supp (w^- - A,)+ C Sfl_i(0) y 1 = 1,2... N 
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and define the functions 

r 1 if Ixl < i? 

X(x) = I and Xi{x) = xix - CO . 

[O if \x\ > R 

Thus we want to find a solution to the problem 

N 

A{W^ + 0) + 5] - A, + </> + i \x'\X x^ = 



in 



with lim|2,|_j.o2 (I'i^) = Oj that is we want to solve the problem 

N 
i=l 

where 

AT 



E■.^AW, + J2{W,-X, + ^u;^\x'\%x^, 



4. A linear theory 

The purpose of this section is to develop a solvability theory for the operator 

N 
i=l 

To this end we introduce the norms 

= sup [ E 1^ - i^\ + 1 ) l</'(^)l , mU = sup IE 1^ - ^^\' + 1 ) IH^)] . 

We want to solve problems of the form L[0] = h with h and having the above 
norms finite. Rather than solving this problem directly, we consider a projected 
problem of the form 

N 3 

L[0] =/i + ^^c,, Z,,x. , (19) 
i=i ]=i 

lim 0(x) = , (20) 
where := dxjWi^ subject to orthogonality conditions 

[ ^Z,,x^dx = Wi = l,2...N, j = 1,2,3. (21) 
Jk3 
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Equation (jl9p involves the coefScients as Lagrange multipliers assoeiated to 
the constraints ((2T|) . If we can solve the equations L[?/;] = h and \-[Yij] = Zij, and 
if we define Cij such that J^^ ^ Zij dx + Cij Jg,, Yij Zij dx = 0, then we observe 
that (/) ^ tp + j Cij Yij solves (fT9)) and satisfies (PT|) . However, for existence, 
we will rather reformulate the question as a constrained variational problem; see 
Equation (gS]) below. 



Lemma 2. Assume that holds. Given h with < +oo, Problem 

(|21[) /las a unique solution (p =: T[ft,] and f/iere exists a positive constant C, 

which is independent of ^ such that, for cu > small enough, 

\m*<C\\h\U,. (22) 

Proof. In order to solve (fTO)) -(pT |) . we first establish ([^^ as an a priori estimate. 
Assume by contradiction the existence of sequences a;„ — ^ 0, ^" satisfying (jl6p 
for w = cj,!, of functions 0„, /i„ and of constants for which 



/ 



||(/)„||* = 1 , lim ||/i„||„ = , 

N 3 

^iiXicfa; = yi,j and L = /t„ + c",- Zjj Xi ■ 

Testing the equation against Zki, we obtain, after an integration by parts, 
0^2/3) |10„||, = / h„ Zu dx + c^, / xfc + ^ |c«.| . 

From these relations we then get c^; = 0(a;,^i/^) + 0(||/i„||„) ^ for all k, I. 

Let us prove that lim„_i.oo ||'?!'ri||L°°(R3) = 0. If not, since ||(/)n||* = 1, we may 
assume that there is an index i and a sufficiently large number i? > for which 

liminf ||(/)„||i^.(B^(5^)) > . 

Using elliptic estimates, and defining tpnix) = 0n(5f +a;), we may assume that ipn 
uniformly converges in sense over compact subsets of M'^ to a bounded, non- 
trivial solution ■(/; of the equation 

^Pd,^w^'xdx==0 Vj = l,2, 3. 

According to [141 Lemma 5], iJj must be a linear combination of the functions 
dxjW'^\ j = 1,2,3. The latter orthogonality conditions yield -0 = 0. This is a 
contradiction and the claim is proven. Finally, let 

N 3 

hri ■■= hn+^Yl ^ij ■ 

Then we have that 

k 
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and hence (/>„, the unique solution in M.'^ of 

^(i>n = hn , hm 0„(x)=O, 

|a:|— ^oo 

satisfies 

\U^)\ < [0{u;'J') + 0{\\hj^.)) ^ . 

Now, since c/)„ — is harmonic in M'^ \Uii3fl(^f ), it tends to zero as |a;| — > oo and 
gets uniformly small on the boundary of this set. By the maximum principle, we 
get the estimate 

\cbni^)\ < {Oiu^'j') + 0{\\hj^.)) ^ on \ U,Bnm . 

This shows that lim„_j.oo ||0n||* = 0, a contradiction with = 1, and (|22l) 

follows. 

Now, for existence issues, we observe that problem (|19I) - (PT|) can be set up in 
variational form in the Hilbert space 

H = |(/)G2?^'2(R^) : / (f>Z,jX^dx = V i = 1, 2 . . . A^, j 1, 2, sj 
endowed with the inner product (0, ^p) = J^., VtJj ■ V0 dx, as 

/ Vcjj-Vtlj dx~ [ p{W^- \i + ^uj^\x'\^Y~\^(t)i} dx+ [ il)hdx = Q 

(23) 

for all t/j ^ H. Since the potential defined by the second term of the above equality 
is compactly supported and h decays sufficiently fast, this equation takes the 
form (p + K[(j)\ = h where K is a compact linear operator of H . The equation for 
h = has just the trivial solution in view of estimate (P^ . Fredholm's alternative 
thus applies to yield existence. This concludes the proof of Lemma [2l □ 

We conclude this section with some considerations on the differentiability of 
the solution with respect to the parameter ^. Let us assume that h ~ 
defines a continuous operator into the space of functions with finite || • |j**-norm. 
We also assume that ^)||** < +oo. Let cj) = 'Pi'iO be the unique solution 

of Problem p^ - (PT|) for that right hand side, with corresponding constants 
Cij{^). Then is diffcrcntiable in ^. Moreover d^cj) can be decomposed as 



d^(j) = ^/j + Y^ d,j Z,j Xj 
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where ijj solves 



N 3 



- Y <9? i^ij X^) + hj Z, 



ij Xi 



1=1 j=l 



lim ip{x) = 

1x1— foo 



t/jZ^jXidx^O V i = 1, 2 . . . N, j = 1, 2, 3 , 



and the constants are chosen so that 77 := X^i^i '^u satisfies 

Xij Zij T] dx ^ - {xij Zij) 4>dx y i = l,2...N, j = 1,2,3 



Lemma 3. With the same notations and conditions as in Lemma\^ we have 

< C(!|M-,0II** + II9?M-,0IU*)- 

5. The projected nonlinear problem 

Next we want to solve a projected version of the nonlinear problem psp us- 
ing Lemma [2] Thus we consider the problem of finding (p with \\4>\\* < +00, 
solution of 

N 3 

l[^]=-E-m+YT. (24) 

i=l ]=l 



lim (/)(a;) = 

|a;|— f+oo 



(25) 



where the coefficients Cy are Lagrange multipliers associated to the orthogonality 
conditions 



.Zy Xidx = Q V ^ = 1, 2 . . .iV, J = 1, 2, 3 . 



(26) 



In other words, wc look for a critical point of the functional J defined by (|13p 
under the constraints ([^5]). 

For this purpose, we first have to measure the error E. We recall that 



N 



1=1 

N 



1=1 



^ p \uH - A, + i ( ^ + i 
1=1 jy:j 



p-1 

+ 
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for some function t taking values in (0, 1). It follows that 

TV N 

from which we deduce the estimate 

|!E||„ <Cc^2/3 

As for the operator N[0], we easily check that for \\(t)\\* < 1, 

N 

I N[0] I < C ^ l^/^r X^ with 7 min{p, 2} , 

1=1 

which implies 

\\m\L<cmz. 

Let T be the linear operator defined in Lemma [2l Equation (|24p can be rewrit- 
ten as 

c^ = m :=-T[E + N [(/.]] . 
Clearly the operator A applies the region 

6= {(/.e L°°(R3) : II0II, < A-a;2/3| 

into itself if the constant K is fixed, large enough. It is straightforward to check 
that N [(j)] satisfies in this region a Lipschitz property of the form 

||N[0i]-N[(/)2]||^^ <K^||01-(/)2||* 

for some positive k^^ such that lim(j_5.o n^^ ~ 0, and hence existence of a unique 
fixed point of A in S immediately follows for cj small enough. We have then 
solved the projected nonlinear problem. 

Since the error E is even with respect to the variable x^, uniqueness of the 
solution of ((24)) - ((26l) implies that this symmetry is also valid for 4> itself, and 
besides, the numbers c^a are automatically all zero. Summarizing, we have proven 
the following result. 

Lemma 4. Assume that ^ = (^i, ^2, • ■ ■S.n) G is given and satisfies ([T5)) . 
Then Problem (j24p - (j26p has a unique solution 0j which depends continuously 
on ^ and io for the \\ \\^,-norm and satisfies H^^H* < CuP'^'^ for some positive C, 
which is independent of u, small enough. Besides, the numbers Ci^ are all equal 
to zero for i = 1, 2 . . . N . 

It is important to mention that (f>^ also defines a continuously differentiable 
operator in its parameter. Indeed, combining its fixed point characterization 
with the implicit function theorem and the result of Lemma [21 we find in fact 
that 

We leave the details to the reader. 

With the complex notation of Section [U let us consider the rotation e*" of an 
angle a around the x'^-axis and let e* " ^ — (e* " ^1 , . . . " ^^r ) . By construction, 
there is a rotational symmetry around the a;'^-axis, which is reflected at the level 
of Problem ([Mll-dMl) as follows. 
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Lemma 5. Consider the solution (f> found in Lemma^ For any a G R and any 

(cc'jCC"^) G C X R, we have that 

The proof is a direct consequence of uniqueness and rotation invariance of the 
equation satisfied by 0^. 



6. The variational reduction 

We consider the functional J defined in (IT51) . Our goal is to find a critical point 
satisfying (^5]). of the form u = + (j)^. We estimate J [W^] by computing first 

/ \VW^\'^ dx = [ IVVwj dx^^ [ |Vw,pdx+V / Vw^ -Vwj dx . 
Jr3 I Jr3 7r3 

The last term of the right hand side can be estimated by 



dx 



- / Aw^ 

JR3 




dx = {wi 

JR3 


JR3 


X^ 




JR3 




'+ 47r|x + e, 


mi rrij 




0(l.4/3) 







where = Jj^j (w^' — '^0+ ~ •^i'^ Next we find that 

JR3 ^ ^. ^ + 



- Xi)^^ dx 



+ {p+l) K-A,)^fVw;j + iw'|a;'|2) + 0(c^4/3) 
^R3 V^.^^ / 



Let us define 



1 f \Vw\^dx-^ ( {w~l)l^'dx 

2 Jr3 p+l Jr3 



Combining the above estimates, we obtain that 

N 



Sm ^^Ar^e, - V-(Ci,...^A.) + 0(c.^/^), (27) 
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where V.|^j(^) = J2i^j stt";™^ I + S'^^ J2iLi ICiP has been introduced in Sec- 
tion [H Here the 0{u}'^/^) term is uniform as w ^ on the set of ^ satisfying the 
constraints (jl6p . This approximation is also uniform in the sense. Indeed, we 
directly check that 

VeJ[M/d = -v«v;;;(C) + o(c.4/3). (28) 

According to ([5]), we have Vl^^^iS,) = uj'^I'^ Vm(C) for C = t^^^^ ^- We get a solution 
of Problem (fTT|l - (fT^ as soon as all constants cij are equal to zero in 

Lemma 6. With the above notations, Cij ~ for all i = 1, 2 . . . j ~ 1, 2, 3 
if and only if is a critical point of the functional ^ i— > A{£^) :~ J [W^ + 4>^]. 

Proof. We have already noticed in Lemma that the numbers c^a are all equal 
to zero. On the other hand, we have that 

/ |Zyfx.dx] + ( J2 cm) 0(^2/3)^ (29) 

From here the assertion of the lemma readily follows, provided that uj is suffi- 
ciently small. □ 

Remark 1. An important observation that follows from the rotation invariance of 
the equation is the following. Assume that the point <^ is such that = (^^i, 0) ^ 
(0, 0) for some ^ G {1, ... TV}. Then if 

%^.yl(0 = for aU fc = 1, . . . , .7 = 1, 2, {k,j) ^ (£, 2) , 

it follows that actually ^ is a critical point of A. Indeed, differentiating in a the 
relation yl(e*"^) = A{S_) we get 

N 

fe=l 

and the result follows. 



7. Proofs of Theorems [TJH 

Let us consider the solution 0^ of ((24|) - ((26)l . i.e. of the problem 

N 3 

lim 0c (a;) = 

|a:|— ^+oo 
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given by Lemma HI We will then get a solution of Problem (TTlTl - lfTS)) . of the 
desired form u = + (f)(_ , inducing the ones for Theorems [1] and [3l if we can 
adjust ^ in such a way that 

Ci] (0=0 for aU i = 1, 2 . . . A^, j = 1, 2, 3 . 

According to Lemma El this is equivalent to finding a critical point of the func- 
tional 

We expand this functional as follows: 

1 

][W^]= S[W^+(ly^]-DS[W^+4>^]-cp^ + - D^S[W^ + {l-t)(t>^]-{(t>^,4>^)dt. 

By definition of 0^ we have that _DJ \W(^ + 0^] • 0^ = 0. On the other hand, using 
Lemma IH we check directly, out of the definition of (j)^ , that 

uniformly on points S,i satisfying constraints (|16p . Hence, from expansion (|27p 
we obtain that 

i=i 

We claim that this expansion also holds in the C"^ sense. Let us first observe 
that 

/ E9^Vy^ dx = V^J [W^] and = Z,j . 

Then, testing equation (j24[) against Zij, we see that 

/ {U[(j)\Zij +L[Zij]) (j) dx ^ ~ I EZ.,j dx + y2cki ZijZuXidx. 

Next we observe that 
\\^[Z^^\l,=O{L0^"^) and / Z,,ZklX^dx = 0{u:^"'') ii (ij) ^ {k,l) . 

By lemma m \\(t>^\\* ~ 0(cl>^/'^), and so we get 

Hence, according to relation (PH)) . we obtain 

(l3jv + 0(w'/')) V4^(0 = V^Jiiy^] + 0{uj^/^) 
where V^J [W^] has been computed in psp . Summarizing, we have found that 

V^AiO = -V^V,'r,(6,-..eiv) + 0(^'/') . 

Therefore, setting ^ = w^/^ C with C = (d, . . . (at) and defining r{() yl(0 on 
03^ := |( g R2Ar . (jjYj) hoidsj^ .^^e t^q^q shown the following result. 
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Proposition 1. With the above notations, we have that 

N 

vr(c) = -vv„,(C) + o{uj^/^) 

uniformly on C satisfying (jl7p . Here the terms O(-) are continuous functions 
of C, defined on 58^. 

7.1. Proof of Theorem]^ If /i > is fixed large enougli, we have that 

inf Vm < inf V,„ . 

Fixing such a we get from Proposition [l]tliat, for all sufficiently small w, 

sup r > sup r . 

so that the functional A has a maximum value somewhere in oj^^^ *B^i , which is 
close to a maximum value of V^. This value is achieved at critical point of A, and 
hence a solution with the desired features exists. The construction is concluded. 

□ 



7.2. Proof of Theorem When (d, . . .(n) is a regular polygon with Q given 
by ^ and all masses are equal, the system is invariant under the rotation 
defined by 

x = {x\x^,x^)^ {{x^ +ix^), x^)^ (e^'^^^ix^ +ix^), x^) ^-.TZnx. 

6R3 gCxR 

We can therefore pass to the quotient with respect to this group of invariance and 
look for solutions u which are invariant under then action oUZn and moreover 
symmetric symmetric with respect to the reflections {x^,x'^,x'^) ^ {x^ ,—x'^ ,x^) 
and (x^,x^,a;^) ^ {x^ ,x^ .—x^). Here we assume that (Ci,...CAr) is contained 
in the plane {x^ = 0} and Ci ~ (r, 0,0). Altogether this amounts to look for 
critical points of the functional 

JiM^l f \Vu\Ux-^f {u-X, + ^Lo'\x'\^Y^\,dx 

where = {x = {x' ,x^) e ^ x E : x' = re^^ s.t. - f < 6* < f } and 

u £ H^{[2) is invariant under the two above reflections and such that Vu -71 = 
on df2 \ {0} X M and xi is the characteristic function of the support of pi ~ 

(u — Ai + Xi in f2i. Here n = n{x) denotes the unit outgoing 

normal vector at x E dQi. With J defined by it is straightforward to see 
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that J [u] = Ji[u] if u is extended to M.^ by assuming that u{TZn x) ~ u{x). 
With these notations, wc find that 

The proof goes as for Theorcm[TJ Because of the symmetry assumptions, cij = 
if j = 2 or 3. Details are left to the reader. □ 



7.3. Proof of Theorem Wc look for a critical point of the functional F of 
Proposition [1] in a neighborhood of a critical point ( of Vm, which is nondegen- 
erate up to rotations. With no loss of generality, we may assume that Ci 7^ 0, 
(i2 = and denote by V„ the restriction of Vm to (R x {0}) x (R2)A^-i g ^. 
Similarly, we denote by f the restriction of T to (R x {0}) x (R^)^-'^. 

By assumption, is a non-degenerate critical point of Vm, i.e. an isolated zero 
of VV„i. Besides, its local degree is non-zero. It follows that on an arbitrarily 
small neighborhood of that point, the degree for is non-zero for all sufficiently 
small Lu. Hence there exists a zero G (R x {0}) x (R^)^"^ of Vr as close 
to C as we wish. From the rotation invariance, it follows that is also a critical 
point of F. The proof of Theorem [3] is concluded. □ 



Appendix A. Facts on Relative Equilibria 

In this Appendix wc have collected some results on the A'^-body problem intro- 
duced in Section [1] which are of interest for the proofs of Theorems [T][31 with a 
list of relevant references. 



Non-degeneracy of relative equilibria in a standard form. Relative equilibria are 
by definition critical points of the function V„i : R^^ — > R defined by 



N , N 



Here we assume that N > 2, and > 0, z = 1, . . . are given parameters. 



Following Smale in we can rewrite this problem as follows. Let us consider 
the (27V — 3)-dimensional manifold 

N 

Sm-={q=iqi,---qN)eM.^^ : ^ (<7,, ^ - (0, 1) , ^ if i ^ j}. 

i=l 

The problem of finding critical points of the functional 

TV 



Um{qi,-.-qN)^^ J2 
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on Sm is equivalent to that of relative equilibria; see for instance [Hj . Let us give 
some details. Let g be a critical point of J7„i on Sm- Then by definition, there 
are Lagrange multipliers A G M and /x G for which 

1 ^ 

~5r2^ \n -n\3 =>^Tnjqj + mjn y j = i,...N . 



First, adding in j the above relations and using that M = Y:.^^ nij > we 
obtain that fj, = 0. Second, taking the scalar product of against qj and then 
adding in j, we easily obtain that UmiQ) = A. From here it follows that the point 
^ = A^/"^ g is a critical point of the functional Vm, hence a relative equilibrium. 
With the reparametrization of given by 



Cia,p,q) = (Cu-.-Cn) = (agi +p,...aqN + p) , {a,p,q) e R x x S', 



the Hessian matrix of Vm at the critical point ( = ({\^^^ ,0,q) found above is 
represented as the block matrix 

\ X-^'DljUmiq), 

where I2 is the 2x2 identity matrix and Ds^^ represents the second covariant 
derivative on Sm- Reciprocally, we check that a critical point 1^ = {Cj)jLi of Vm 

necessarily satisfies J2^=i ^''■j (j = 0- Defining q = (^X^jLi'^j 10 P) C we 
readily check that g is a critical point of Um in Sm- 

Any rotation e' " g of a critical point q of Um on Sm is also a critical point. 
We say that two such critical points are equivalent in 5"^. Let us denote by Sm 
the quotient manifold of 5,„ by this equivalence relation. On Sm, critical points 
of the potential Um yield critical points of Um on Sm and hence equivalence 
classes of critical points e* " C for Vm using the reparametrization. 

A critical point q of Um on Sm is said to be non-degenerate if the second 
variation of Um at q is non-singular. Let us assume that qt ^ 0, with either 
= 1, or £ = 2 if gi = 0. Then there is a unique representative g of this class of 
equivalence for which qi2 = 0. It is a routine verification to check that g is then 
a critical point of Um on the (2A^ — 4)-dimensional manifold 

S,„ ;= {g e Sm ■ qi ^0 as above , g^2 = 0} . 

Moreover, the second derivative of [/,„ on 5,„ at g is non-degenerate if and only if 
Um{q) is non-singular. Because of the expression of D'^VmiO, '^^ see that C, 
is non-degenerate as a critical point of V,„ on the space of C G R^^ with qg2 = 0, 
which is the notion of non- degeneracy up to rotations of a relative equilibrium 
that we have used in this paper. Finally we define the index of a non-degenerate 
relative equilibrium C, as the number of negative eigenvalues of Um{q)- 
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Some results on classification of relative equilibria. For simplicity, wc will assume 
that masses are all different: for any i, j = 1, . . . N, if nii = m^, then i = j. This 
is the generic case. 

The cases iV = 2, 3 are well known; sec for instance For N ~ 2, the only 
class of critical points is such that 

iCi - C2I = (If)^^^ and miCi+TO2C2 = with = mi + ma . 

For TV = 3, there are two types of solutions, the Lagrange and the Euler solutions. 
The Lagrange solutions are such that their center of mass is fixed at the origin, 
the masses are located at the vertices of an equilateral triangle, and the distance 
between each point is (Af/(47r))^/'^ with M = mi + m2 + ma. They give rise 
to two classes of solutions corresponding to the two orientations of the triangle 
when labeled by the masses. The Euler solutions are made of aligned points and 
provide three classes of critical points, one for each ordering of the masses on 
the line. 

In the case > 4, the classes of solutions for which all points are coUinear 
still exist (see [55]) and are known as the Moulton solutions. But the config- 
uration of relative equilibria where all particles are located at the vertices of 
a regular A^-polygon exists if and only if all masses are equal; see [26ll43ll35l 
[T^H5] . Various classification results which have been obtained by Palmore are 
summarized below. 

Theorem 4 ( [30l[3n[32l[33l[34] 1 . We have the following multiplicity results. 

(a) For N > 3, the index of a relative equilibrium is always greater or equal than 
N ~- 2. This bound is achieved by Moulton' s solutions. 

(b) For N > 2,, there are at least Hi{N) := {^){N ~ 1 ~ i){N - 2) \ distinct 
relative equilibria in Sm of index 2N — 4 — i if Um *s o, Morse function. As a 
consequence, there are at least 

J2 = [2^-i(7V - 2) + 1] (iV - 2) ! 

i=0 

distinct relative equilibria in Sm if Um is a Morse function. 

(c) For every N > 3 and for almost all masses m S M.^ , Um is a Morse function. 

(d) There are only finitely many classes of relative equilibria for every N > 3 and 
for almost all masses m = {'mi)iLi G IR+ • 
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